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We present a phenomenological investigation of non-Gaussian effects that could be seen on CMB 
temperature maps. Explicit expressions for the temperature correlation functions are given for 
different types of primordial mode couplings. We argue that a simplified description of the radial 
transfer function for the temperature anisotropies allows to get insights into the general properties 
of the bi and tri-spectra. The accuracy of these results is explored together with the use of the small 
scale approximation to get explicit expressions of high order spectra. The bi-spectrum is found to 
have alternate signs for the successive acoustic peaks. Sign patterns for the trispectra are more 
complicated and depend specifically on the type of metric couplings. Local primordial couplings are 
\Q • found to give patterns that are different from those expected from weak lensing effects. 
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I. INTRODUCTION 



The standard inflationary models [| have been very successful in explaining the basic features of the CMB observa- 
tions: near scale invariant power spectrum of, to a good approximation, primordial Gaussian adiabatic perturbations. 
These are undoubtedly important results that give support to the inflationary scheme. From another point of view it 
is an uncomfortable situation since there is practically no possibility with the current data sets to distinguish different 
J> models of inflations. The search of observational signatures beyond predictions from generic inflation is therefore 
crucial for getting insights into the nature of the inflaton field or more generally into physics of the early Universe. 
For instance trans-Planckian effects 0, presence of isocurvature modes 0, 0, and non-Gaussian effects 6] could 
betray some aspects of the inflationary physics. It is now well established that in standard inflation, e.g. single field 
inflation with slow roll conditions, a minimal amount of non-Gaussianities is expected to be induced along the cosmic 
evolution [fjj. They simply come from the couplings of modes contained in the Einstein equations, that are intrinsi- 
cally nonlinear in the fields. The more famous of those are the lensing effects on the CMB sky They indeed 
can be seen as the effects of couplings between the gravitational potential on the last-scattering surface and those 
present on the line-of-sights. Of course other local couplings exist that take place before or during the recombination 
period. For instance Bartolo and collaborators in [(| derive the amplitude and nature of the superhorizon couplings 
of the gravitational potentials in terms of their expected bi-spectrum. Its shape betrays the nature and shape of the 
i-j ■ quadratic couplings in the Einstein equations that generically are expected to determine the types of couplings for 
^ \ the modes that reenter the horizon. 

For one-field inflation couplings at horizon crossing are rather generic and induce typically 10~ J effects in the 
metric perturbations. Effects of such amplitude are expected to be marginally visible. They are however expected 
to be altered by second order effects taking place at low redshift after horizon crossing. Except for the Sachs Wolfe 
plateau or the lens effects, predictions of the observable intrinsic non-Gaussian effects are not known. They require 
the computation of the physics of recombination up to second order. An enterprise still to be done. 

There are however cases where significant mode couplings could survive the inflationary period. This is the case for 
some flavors of the curvaton models This could also be the case in other models of multi-field inflation where a 
transfer of modes, from isocurvature to adiabatic, is possible during or at the end of the inflationary period as described 
in 0, HTl IT^ | . In this case no isocurvature modes are expected to survive, contrary to the case of the curvaton model, 
and the models predict only adiabatic fluctuations the spectrum of which can be arbitrarily close to scale invariancc. 
However nothing prevents the initial isocurvature modes from developping significant non-Gaussian features which 
may be transfered into those of the adiabatic fluctuations. Such a transfer mechanism has been described in details 
in [l3| where the resulting high order correlation properties are explicitly given. 

In principle, the properties of the metric fluctuations entirely determine those of the temperature or polarization 
CMB sky. However, the transcription of the properties of metric to properties in the observed temperature field 
reveals involved for at least two reasons. First, it is to be noted that the transfer physics is naturally best encoded in 
Fourier space where all modes evolve independently from one-another whereas the nonlinear couplings of Eq. J2| is 



'Electronic address: Tristan.Brunier@cea.fr 
^Electronic address: Francis.Bernardeau@cea.fr 



2 



local in real space. The second reason is that projection effects should be taken into account so that the modes that 
are observed correspond to a collection of Fourier modes. The then simple functional forms (such those of Eqs. Q 
and (JHJ for instance) that are rather generically expected for the metric are then greatly altered. For instance the 
sign of the three- and four-point correlation function may depend on scale. The aim of the investigations pursued 
in this paper is to try to uncover simple prescriptions for the shape, e.g. angular dependence, of those quantities. 
Comparison with the trispectra induced by weak lensing effects will also be presented. 

The paper is divided as follows. In the second section we detail the couplings we assume for the primordial metric 
fluctuations. In the third section we present the basic quantities that are required to describe the way the potential 
fluctuations are transfered into the temperature fluctuations. We then take advantage of these results to present 
the theoretical shapes of the bi- and tri-spectra of the temperature field in case for the models of potential high- 
order correlation functions described in the previous paragraph. The following section is devoted to an attempt to 
obtain a phcnomcnological description of the transfer function that gives insights into the angular dependence of these 
correlation functions. The last section is devoted to computation in the small angle approximation. Calculations are 
obviously much more straightforward in this limit which indeed corresponds to a regime where most of the observations 
can be made. 



II. A MODEL OF PRIMORDIAL NON-GAUSSIAN METRIC PERTURBATIONS 



In the analysis pursued in this paper we assume that the primordial metric fluctuations are those expected in rather 
generic models of multiple-field inflation, more specifically along the descriptions presented in [Tol ITU Il4|. In this 
family of models, the surviving couplings in the metric, are expected to be, to a good approximation, equivalent to 
those induced by the superposition of two stochastically independent fields, a Gaussian one and one obtained by a 
non-linear transform of a Gaussian field with the same spectrum. In other words the local potential would read, 

= cos a $i(a?) + sin a J r [$ 2 (x)] (1) 

where ratio between the initial adiabatic and isocurvature fluctuations is described by the mixing angle a and where 
the function T takes into account the self coupling of the field that gave rise to this part of the metric fluctuations. The 
function T obviously depends on the details of the inflationary model in particular on the self interaction potential 
along the transverse directions. As argued in |10| a natural choice is a quartic potential. In such a case, the function 
T is characterized by two quantities. One is the amplitude of the coupling constant times N e , the number of e- foldings 
between horizon crossing of the observable modes and their transfer from isocurvature to adiabatic direction, the other 
is related to the field value in the transverse direction $ at our current Hubble scale. The latter actually corresponds 
to a finite volume effect ^4|. In this framework the function T then reads |32l| 

$ + $ 

^($) = , (2) 



'l-|i3($ + *) a /3 

where /13 is related to the self coupling of the field. For a coupling in A0 4 /4 we have 

^3 ~ -AiV e /C 2 2 (3) 

where C2 is the amplitude of the metric fluctuation at super-Hubble scale, C2 ~ 10~ 5 . One important feature of this 
description is that the non linear transform of the field is local in real space. This description has the advantage of 
providing a full description of the model. This equation should however be used with care. It gives a good account 
of the mode couplings only when the effective coupling constant XN e is small. It is therefore preferable to consider 
the bi- and tri-spectra of the potential field. So let us define the spectrum P(k), bi- and tri-spectra, B(k\,k2) and 
T(fci, lt2, ^3), of the potential field the following way, 

= <5 D irac(£l + k 2 ) P(fci) (4) 
(Sjfj = ^Dirac(fcl + kl + h) B(h,k 2) k 3 ) (5) 

. . . % 4 ) c = <5 D irac(fcl + . . . + ki) T(kl,k2, ki, fa), (6) 

where the underscript c stands for the connected part of the corresponding ensemble average. As $1 and $ 2 have the 
same power spectrum, if the nonlinear couplings of $2 are small then the power spectrum of the metric fluctuations 
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is left unchanged (to corrections that are quadratic in the coupling parameter). Then for the form presented in Eq. 
PJl. for a small coupling parameter fj, 3 we have, 

B(k u k 2 ,k 3 ) = v 2 [P(h)P(k 2 ) + perm.] (7) 

with v 2 = sin 3 a /i3$, $ being a priori of the order of H being the value of the Hubble constant at horizon crossing, 
and 

T(k u k 2) k 3 ) = ^ tar [P(fc 1 )P(fc 2 )F(fc 3 ) + perm.] + iJ™ [p(fc 1 )P(fc 2 )P(|fc 2 + fc 3 |) + perm.] (8) 

with f| tar = sin 4 a /i 3 and ^ 3 lnc = sin 4 a fJ,^ 2 . 

Note that if we are interested in the bispectrum, the parameter v 2 is identical to the parameter /jvl usually used 
to describe the nonlinear transform of the metric fluctuations |l6j . It is to be noted however that whereas /jvl is 
expected of the order unity in single field inflation, v 2 can reach much larger values in multiple-field inflation. This is 
one justification of the investigations presented in this paper. 

We can see that the tri-spectrum contains two terms of different geometrical shapes. The relative importance of 
these two terms depends in particular on the value of <&. This ratio is naturally of order unity, but otherwise totally 
unpredictable |33j ] . In the following we therefore consider the two cases and examine the consequences of both terms 
on the temperature high order correlation functions. 



III. EXPRESSIONS OF THE CORRELATION FUNCTIONS 



In this section we set the basic relations that relate the primordial metric fluctuations to the temperature fluctua- 
tions. 



A. Radial transfer function 



The observed temperature anisotropies are decomposed in a sum of spherical harmonics with coefficients ai m . The 
inverse relation gives the expression of those coefficients as a function of the temperature fluctuation ^(7) in the 
direction 7 on the sky, 

a lm = J d 2 7^(7)^„(7)- (9) 

When the linear theory is applied to the metric and density fluctuations the coefficient at m are linearly related to the 
primordial metric fluctuations. In the following we will only take into account the existence of adiabatic fluctuations. 
Then one can write, 

a lm = M-i) 1 J ^^T e (kMk)Y; m (k), (10) 

where $(fc) is the three-dimensional Fourier transformed of the gravitational potential, T^(fc) is the photon transfer 
function and Y^ m (k) are the spherical harmonics in the direction given by the unit vector k = k/k. 

The functions Tt(k) encode all the micro-physics that takes place after horizon crossing until the photons reach the 
observer. There is one well known limit case which corresponds to large-angular scales for vanishing curvature and 
cosmological constant. In this limit indeed, the observed local temperature fluctuation is simply one third of the local 
potential. This is the so-called Sachs- Wolfe effect In this case we simply have Tg(k) = — I/3j^(fcr») neglecting 
the late and early Integrated Sachs Wolfe effects. Although the validity regime of this form is limited, it is worth 
investigating since it corresponds to a case where only projections effects have to be taken into account. 

One aim of this paper is to take advantage of the expression pOI) to relate as accurately as possible the statistical 
properties of the ai m to those of the potential field. To do so we found fruitful to introduce the radial transfer function. 
Let us first decompose the gravitational potential intercepted on a sphere of radius r as 



$(r 7 ) = ^^ m (r)Y, m ( 7 ). 



(11) 
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Note that conversely 

*ft»(r) = J d 2 7$(r7)r/ m (7). (12) 

In Fourier space, we then get 

*® = J ^me rk - S , (13) 

= 7^71 E / r " dr d ^ (14) 

Using the decomposition of plane waves on a sphere 

e ** = 4?r ^ i L j L (kr)Y LM (k)Y£ M (% (15) 

LM 

with x = ry, one obtains 

$(fc) = 4^ J] i L | -T^y^r 2 1 d^^^yL^^y^W^CO^^), (16) 



L.l.M.m 



= y^E^/ drr 2 ^(fcr)y m (fc)$ £m (r 



(17) 



Conversely 



7T 

Hence Eq. Hl()|) can be recast in 



= (-t)V- / d 3 fc$(% £ (fcr)F/ m (fc). (18) 



ai m = - i drR e {r)<P em (r) 7 (19) 
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where the radial transfer function Ri(r) is obtained through a radial Fourier transform of the initial Fourier space 
transfer functions, 

R e (r) = J dkk 2 r 2 T e {k) Je (kr). (20) 

Obviously the micro-physics of the recombination can equally be described by the functions Tg(k) or by the functions 
Re{r). However the latter are more directly related to the physics of recombination, which takes place over a short 
range of radius, and one can therefore try to find approximate forms for their r dependence. The Sachs Wolfe regime 
corresponds to a limit case for which, 

flf W M=-J<W ac (r-r*). (21) 
6 

Before we propose approximate forms for Ri(r) let us first detail how the statistical quantities we are interested in 
are related to this function. 

B. Power spectrum 

We first want to express the ensemble average of products of ag m ag, m , as a function of the potential power spectrum 
and the radial transfer function. Using Eq. f2H|l we naturally get 

(a em a* e , m> ) = feW Q J drr 2 J dr' r' 2 J dkk 2 P(k)Mkr)jr(kr') 

dkxklT^uikxr) [ dk 2 k%T e ,(k 2 )j e ,(k 2 r'), (22) 



where the integration over r' may be performed if one takes advantage of 

da; x 2 jt(kx)je(k'x) = -^S Dirac (k - k'). (23) 

Then temperature anisotropy power spectrum Ct, defined as, 

{ae m a* t , m ,) = (—1)™ {at m at'-m') = &ll'&mm'Ci , (24) 

reads 

°t = [I) J <irR e (r)Ur) , (25) 

where we define 

&(r) = J dkk 2 P(k)T e (k) Je (kr). (26) 
We will see now that the functional relation (|2*51) can be generalized to higher order correlation functions. 

C. Bispectrum 

1. Full sky expression 

We are interested in the ensemble average of the product of three ae m that should not vanish when the potential 
field exhibits non-Gaussian properties. This correlation function formally reads, 

2 X3 



(ae imi ae 2m2 ae ama } = I -J j dn J dr 2 J dr 3 R ei {n)Ri 2 (r 2 )P^ 3 (r 3 )(^ imi {ri)$e 2m2 (r 2 )$e 3 m 3 (r 3 )) ■ (27) 
The ensemble average that appears in the r.h.s. of this equation can be related to the potential bispectrum, 

($Wri)$ &ro2 (r 2 )^ 3m3 (r 3 )) - (-*)*+*•+*» (^j ' J d% d% d% <*(fci)$(fc 2 )*(fc 3 )) 

3iAkiri)3i 2 {k 2 r 2 ) 3l3 {hr 5 ) Y^k^Y^k^Y^k) . (28) 
The Dirac distribution that appears in JJJ) can be rewritten as a Fourier transform 

ferae + fc 2 + h) = J , (29) 

and then expanded into spherical Bessel functions with Eq. (|15fl 

/j3 -» 
— ^ i L ^ L * +L *j Ll (k l x)j L2 {k 2 x) jL3 {kzx) 

^ ^ L 1 ,L 2 ,L- 3 ,MiM2,Mi 

Y LlMl (h)Y L2 M 2 (k)Y L3M3 (k)Y£ lMl {l)Yt 2 M 2 (l)YZ 3 M 3 (7), (30) 
where x = xj. Using the form (JJJ for the expression of the bispectrum and inserting Eq. (|3L)[) into l|28|l . we get 

($Wn)^ 2m2 (r 2 )^ 3m3 (r 3 )>c = 8. 2 ^j 7 '^(X^i^^jyim&fM 



dfc 2 klj t2 {k 2 r 2 )j t2 {k 2 x) J dk 3 k 2 j e3 (k 3 r 3 )j e3 (k 3 x) 
[P(fc 1 )P(fc 2 )+P(fc 2 )P(fc 3 ) + J P(fc3)P(fci)] , (31) 

where the Gaunt integral is defined by 

(3) ^XiT = / d2 7 ^.(fc^W ■ (32) 
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We integrate over the momenta fcj which do not appear as an argument P(ki) in Eq. (|31|l and then perform the 
integrations over two of the radial variables. Eq. (|27|) then becomes 

/r,\ 3/2 

/ \ oil (3) r>mim2m3 

[ai imi ai 2m2 a£ 3m3 ) c — 8^2 I — I * tibia 

J dr [R tl (r)& (r)& 3 (r) + i?, 2 (r)& 3 (r)& (r) + i?, 3 (r)& (r)& (r)] . (33) 

The bispectrum is expressed in terms of the fundamental functions Re and £e that appeared in the computation of 
the power spectrum. We shall see a diagrammatic interpretation of each of these terms in section Till El For instance, 
the third term of Eq. Ij33(l is diagrammatically represented in Fig. ^ 

The expression (|33|l also involves geometrical factors encoded by the Gaunt integrals. In the following, we introduce 
an estimator which should allow to define a reduced quantity. Wc also present the correspondence with the small angle 
approximated bispectrum. 



2. Estimator 



Following [l6L ITt| an unbiased estimator of the angular averaged bispectrum may be chosen to be 

lx 4 l: 
mi ?TJ2 m 

where the Wigner 3-j symbol is related to the Gaunt integral by 



«i«2<;3 — 

mi ,7712 ,7773 



51 ( m \ ml m 3 ) fl ''^ ffi ^» fl( « ' ( 34 ) 



Using the equality 



(3)^777 7 777 27.3 _ / (2ll + 1) (2l 2 + 1) (2l 3 + 1 ) ( l X l 2 £ 3 \fh 4 4 , 

y ti**t* _ V 4tt I mi m 2 m 3 ) I ' ' 



y- / 4 4 4 \( 4 4 t' a \ = ^Wr^ 
\ mi m 2 m 3 J I mi m 2 mg J 24 + 1 ' 



(24 + 1) (24 + 1) (24 + 1) fit 4 4 



we get 



4tt V 

with the reduced bispectrum defined as 

■ 2 x3/2 



Ws = 8^ 2 - dr [R tl (r)& 2 M6 3 M + ^ 2 (r)& 3 (r)& M + ^ 3 (r)& (r)& 2 (r)] ■ (38) 



The reduced bispectrum reveals convenient to describe the non-Gaussian part of the signal as it does not include the 
overall geometrical factors. 

It is also convenient to introduce a normalized bispectrum bi u i it i a defined by 

I b eie 2 e 3 r „ Q \ 

otiMs - 7^~~F< — ~Tn~~n — T n r ' ^ 

W1W2 + W2W3 + 

which reads in terms of the functions Ri(r) and £i(r) 

o \ 3/2 

2 \ ^2 



&Wa = 8 - r r r r / dr ^ ( r )^ (0 + Rt 3 (r)& 3 (r) + ft 3 (r)& (r)] 

k tt j <^e 1 L^e 2 + Ce 2 Oe 3 + o< 3 0£ 1 - 



(40) 
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FIG. 1: Diagrammatic representation of the third term in the expression 13311 of the bispectrum. Each line between a non- vertex 
point (£i,irii) and a vertex (r, £3,7713) represents the term ^(r). The vertex value is G™^™^ 3 Ri?,{ r )- The radial vertex 

position should be integrated over and the whole graph should be multiplied by 8 (-|) 3 ^ 2 . 



3. Small angle approximation 

The reduced bispectrum encodes all the physical processes that lead to a non vanishing bispectrum. On the other 
hand, the Gaunt integral that appears in Ea. l|38[) only carries an overall geometrical dependence and ensures that the 
momenta ^1,^2 and £3 satisfy the triangular inequalities. 

It appears that this overall geometrical dependence translates into a simple momentum conservation in the flat sky 
approximation (see Appendix [UJ . Defining the quantity a{£) as in Appendix O the bispectrum reads in the small 
angle approximation 

(a(/ 1 ) a (f 2 ) a (f 3 ))c - -U^^Dirac^l + h + £3) (41) 

where 6^ 2 ^ 3 is the reduced bispectrum defined in Eq. i|38H . The Dirac function imposes that the vectors £i,£2 and 
£ 3 form a triangle whose lengths respectively match with £i,£ 2 and £3. 
We now turn to the study of the tri-spectrum. 



D. Trispectrum 

1. Full sky expression 



2^ 4 



Generally the trispectrum reads 

J j dri J dr 2 J dr 3 J driR£ 1 (ri)R£ 2 (r 2 )Rt 3 (r3)Ru( r '0 

($i imi (n)^ 2 m 2 (r2)*< s m, (rs)^^ (r 4 )> . (42) 

From the expression (JHJ), we will have 2 contributions. One star contribution, {ae imi ■ ■ ■<^ 4 m 4 ) star due to the first 
term and one line contribution, (o^ imi . . . ag inii ) hnc , due to the second one. These denominations will actually become 
much clearer in section Till El where a diagrammatic representation of those terms is developed. For instance, only the 
diagrams whose shape is drawn in Fig. [21 contribute to the line trispectrum, whereas the star trispectrum is made of 
terms whose shape is drawn in Fig. |21 

For the line contribution, the connected four-point function of the potential is given by 

. . . $(fc 4 ))c ino = 4 inc / d 3 fc 5 P{h) [P(ki) + P(k 2 )} [P(k 3 ) + P(k 4 )} <5 Dirac (fci + k 2 - fc 5 ) <5 Dirac (fcs + k 4 + fc 5 ) 
+ (2 <-> 3) + (2 «-> 4) , (43) 
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which entails using Eq. i|30|) . 



2 T 2 
X 2 



«,«P(hW ll (h n )fc(h«)/«it,«fc(t«)fc(»«) 

dfc 3 k 3 P(fa)jt 3 (for 3 )ji 3 (k 3 x 2 ) J dfa kjje 4 (far 4)3^ (fax 2 ) 

dfc 5 k 2 5 P(k 5 )j L {k 5Xl )jL{k 5 X2) + (1 <-> 2) + (3 4) + f J ~ J 

+ (2 <-> 3) + (2 ^ 4) . (44) 
Integrations over fc 2 and fc 4 in the first term impose x\ = r 2 and x 2 = r 4 such as 

(<f£ imi (r 1 )$, 2m2 (r 2 )^ 3m3 (r3)^ 4m4 (r 4 ))i i " c = ^™2 8 £ (-1)^^^^ / d*i tflWi* (fei^ for*) 

L,M ^ 
.2 r>/ 1. "v • /;„ „ \ „• fi. „ \ / ji. 7.2 i 



dfc 3 k 3 P(k 3 )jt 3 (k 3 r 3 )j i3 (k 3 r 4 ) J dk 5 k 5 P(k 5 )j L (k 5 r 2 )jL(hr4,) 

+ 11 other terms. . (45) 

Using the definition of the radial transfer function H2U[) . integrations over two of the radial variables may be 
performed within the expression of the trispectrum which becomes 



/_ _ _ _ \linc g2 line f* 1 \ I i \ M s1m1m.2MsTm.3m4— M 

{ae 1 m 1 ae 2 m 2 a,£ 3 m 3 ae i m i ) c - s v z |— J ^(-J-J ^ 3 ^ 4 z, 

^ ' L.M 

dr J dr'R ei (r)Z i2 (r)t L (r,r')R i3 (r')& 4 (r') + (1 «-> 2) + (3 ~ 4) + ( 



1 <-> 2 
3^4 

+ (2 «-> 3) + (2 «-» 4) , (46) 

where 

C L (r,r') = f dkk 2 P(k)j L (kr)j L (kr') . (47) 



Contrarily to the previous correlation functions, the expression of the trispectrum involves a third function which 
only depends on the primordial power spectrum. We will see in section IIII El that this term has a diagrammatic 
interpretation. For instance, the second term in Eq. (|46|l is diagrammatically represented in Fig. 
For the star contribution we have, 

(^ imi (n)^ 2m2 (r 2 )^ 3ro 3(r3)^ 4m4 (r4)) 8 c tar = f ^Yl^^Y^^Y^^^,) J <^ X 2 

dh kfj £l (far^j^ (fax) J dk 2 k\u 2 [k 2 r 2 )j l2 {k 2 x) 

dfa klji 3 {k 3 r 3 )ji 3 (k 3 x) J dfa k\u 4 (far 4) fa (fax) 
[P(fa)P(k 2 )P(k 3 ) + 3 other terms] . (48) 
The first term should be computed by integrating first over fa, which imposes x — r 4 . This gives 

^^mAri^^mAr^^mAr^^ir^ = 2 6 ^ tar J d 2 ^, (j)Y e * 2m2 (j)Y e * ma (WlmM 

dki klP(k 1 )je 1 (k 1 r 1 )je 1 (k 1 r 4 ) J dk 2 k 2 P(k 2 )j i2 (k 2 r 2 )je 2 (k 2 r 4 ) 

dk 3 kjP(k 3 )j£ 3 (k 3 r 3 )ji 3 (k 3 r 4 ) + 3 other terms . (49) 
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The angular integral may be expressed as a function of Gaunt integrals 



L.M 

Using the same expansion as in the previous subsection, we get 
{a limi a l2m2 a l3mz a lim ^ = ^ tal - £ {3) QTXl M (^"^""(-l)" f dr [R H (r)& a (r)&, (r)& (r) + 3 other terms] 

^ L,M ^ 

(51) 

One of the terms in Eq. (|51|l is diagrammatically represented in Fig. |3J 

The expressions 1461 and 15 II are in perfect agreement with |l9l | in the case of the trispectrum. Indeed, this formalism 
can easily be extended to higher order terms. Any N-point correlation function may be expressed in terms of the three 
functions Ri, ^ and (,l using some basic rules that are detailed in section Till El and associated with a diagrammatic 
description. 

Before going further, we stress how to decouple the overall geometrical dependence given in terms of Gaunt integals 
from the primordial signal. To do so, we build an estimator for the full-sky trispectrum and study its limit in the 
small angle approximation. 

2. Estimator 

Following [Til , Hsj , we write the trispectrum in a rotational invariant form 

{a iimi a i2m2 a t , m3 a Umi ) c = ^(-1) M m 2 m){ m 3 m\ —M ) ' (52) 

L,M x ' x ' 

Similarly to the bispectrum, an estimator for the connected part of the angular averaged trispectrum may be chosen 
to be 



TfltHL) = (2L + E^ 1 ^ ( t rt M ) ( t ml —M ) " G Zt( L )> ( 53 ) 

all m M ^ ' ^ ' 

where G^ 2 (L) is the estimator of the unconnected terms and is defined in such a way that T^ 2 (L) vanishes for a 
Gaussian field (see 0] for a discussion though the observable is not explicitly written in there), 



oft (l) = (2i+i)V Y(-i) M ( £l i2 L)( h iA l m ) 

' v ; ^— i ^ x ' \ mi mo M I \ 7713 m4 — M I 

all m M ' - / \ 





^2 




V mi 


m 2 


£ ( 



r^mi— m2^m 3 -m4^i^2^3-«4 / , ( — 1)™ + "' afim' a £i-m'af 3 m"a£ 3 -m" 



m' ,m'' 



[(24 + l)(2^ 3 + l + 2^ 3 ) 

+ (2^3)7(2^4) L. (54) 



In the star configuration, we define the reduced averaged trispectrum star i^^(L) as 



«~T*£{L) = V(24 + 1)(24 + 1)(24 + 1)(24 + 1) ( o 1 ^ ) ( ) ^'tt^ ■ ( 55 ) 

so that it takes the L-independent form 
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The case of the line configuration is much involved since it implies three kinds of terms with different geometrical 
configurations encoded by the Gaunt integrals. The line trispectrum takes the form 



1 1 10 



L,M 

with 



e« h l 



= El"')" 1 1 M 1 1 -M WW + P » 3) + (2„ 4) , (57) 



f&'.'TO = ^V(2<, + l)(2f s + l)(24+lX24 + l)(^ f 5 ) ( S ''S ) """<«^) - 0») 
and the reduced trispectrum lme ^ 1 ^(L) is given by 

line t^ 4 (£) = 8 2 ^" c (|) 2 { / dr / dr'i?, 1 (r)6 2 (r)a(r,r') J R fe (^')^(^) + (1 - 2) + (3 - 4) 
Using the link between the Wigner-6j and Wigner- 3j symbols (see EUllll^) 




the quantity T^ 1 ^ (L) can be written as 

*~t;;,':<i> = p;,t(i) + (2i+i)i;(-i)''«-{^ \ ^^tm 



+< 2i + i)D-i) i+i " {t t L L «y, 



■1 *4 

3^2 



(L") . (61) 



The Wigner-6j symbol | ^ d / j" re P resen ^ s a quadrilateral with sides (a, 6, c, d) whose diagonal form the triangles 

(a, d, /), (b, c, /), (c, d, e), (a, 6, e) and vanishes if any of the related triangular inequalities are not fulfilled (see p7jl- 

Note that contrary to the star trispectrum, the line trispectrum is generically singular when L = 0. This is possible 
if l\ — li and £3 = £4 (and symmetric cases) and is due to the expression (|59H which is then logarithmically divergent 
for a scale invariant power spectrum. This apparent divergence however does not appear in observable quantities such 
that in Eq. (|53[1 . Indeed it can be easily checked that terms involving Pf^(L — 0) and coming from the two terms 
of i|53|) exactly cancel each other. Super-Hubble effects remain then unobservable. 

Finally, we find convenient to introduce the normalized trispectra (L) defined as 



H t i 2 ^1%U (L) 



{c tl + c l2 )c L {c l3 + c u y 10 ' 



3. Small angle approximation 

Using the notations of Appendix IO we can define the quantities a(£) for large enough multipoles so that the line 
and star trispectra take the forms 

(a(4)a(^ 2 )a(4)a(4)) 8 c t ar = W (h + h + £3 + h) staT t% \\ (63) 

and 

(a(fi)a(/2)a(4)a(f 4 ))c Iie = &ir ac (4 +h + £3 + 4) J d 2 L [5 Dirac (/i + £ 2 - £) line ^(L) 



+folrac(*l + £3 - Lf n %lll(L) + fohac^l + h ~ Lf**1%£ (L)^64) 

The Dirac functions ensures that the multipoles £%,£%, £3, £4 form a quadrilateral. The star-trispectrum does not 
depend on the shape of the quadrangle since it is enterely determined by the side lengths. Yet the line-trispectrum 
not only depends on the side lengths but also on the diagonals L12 = £\ + £2, £13 = £\ + £3 and Z14 = £\ + £4 of the 
quadrilateral formed by the four multipoles. 



11 




FIG. 2: Diagrammatic representation of the "line" contribution to the trispectrum. Each line between a non-vertex point 
(£i,mi) and a vertex (r,£,m) represents the term ^(r). Each line between the vertices (r, £2, "12) and [r 1 ,£3,1713) represents 
the term Ct(r, r'). (Here L is the total angular momentum from l\ and £2 or equivalently from £3 and £4). The vertex values 
are i/ 2 (3) QT^l M R e 2 ( r ) and v ^ 3) GT\Tl~ M - R « 3 ( r ')- Tne radial positions of the vertices should be integrated over and the 
internal momenta (L,M) should be summed over. The whole graph should be multiplied by 8 2 (§) 2 - 




FIG. 3: Diagrammatic representation of the "star" contribution to the trispectrum. Each single line between (£i,m,i) and the 
vertex (r,£,m) represents the term ^(r). The vertex value at the point r is given by QT^i^i^il ™ l -^^M- The radial 
position of the vertex should be integrated over. The whole graph should be multiplied by 2l. 



E. Diagrammatic representations 



The formal results obtained in the previous sections show that at least some aspects of the potential correlation 
functions are preserved. That furthermore suggests general rules for the constructions of the N-point correlation 
functions for the spherical harmonics coefficients. Let us start with basics rules for the potential correlation function. 
Using a usual diagrammatic representation, assume that the following rules hold for any potential iV-point function 
) c for tree order or loop graphs: 

• each vertex connecting rij lines contributes to a factor v ni -i and a Dirac distribution of the sum over all the 
connected momenta, 

• each external line carrying the momentum fcj represents P(ki). For pi external lines connected to a vertex, the 
contribution is a sum over configurations of the product of pi — 1 potential power spectra P(ki), 

• each link between two vertices i and j should be formed by internal line. Each internal line stands for P{qi) and 
imposes a momentum conservation through ^Dirac^ + Qj), where qt (respectively qj) is the internal momentum 
out of the vertex i (respectively j). 

To compute the induced iV-point temperature correlation function, we use the basic relation 

/ 2 \ 3N/2 f 
(a eimi ...a tNmN ) = f-J (-i) il+ - +iN dri...dr N R ei (ri)...R tN (r N ) 
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jt 1 (k 1 ri)...j tir (k N r N ) Y e * imi (k 1 )...Y t * NmN (k N ) . (65) 

In the expression of the potential ./V-point function, the momentum conservation at the vertex i with rii—pi external 
lines and pi internal lines contributes as 

^2 i Li+-+L n . (47r) f d s x 4 j L i {k-[Xi)...j L i^ ^ (fc n< _ P4 x i )j L i^_ () _ + i (qiXi)...j L i ni {q Vi Xi) 

L{,...,L' n .,Mi,...,M^. 

YL*M'(ki)-Y L i M i (k ni - Pi )Y L i M i (qi)..Y L i Mi {q Pi ) 

1 1 n i -p i n i —p i n i -p i + l n$ — p£ + l n i i 

Y* tMi (%)...Yl kMii m , (66) 

where k (q) stands for an external (internal) momentum. 

Integration over the angles will contribute with the overall geometrical factor (ni 'Q Li L ? . Integration over the 

angles fe, in Eq. I|65(l will impose = £j for all external momenta. Integrations over all the external momenta that do 
not appear as P(ki) in the expansion Eq. lj6*5|l give it/2 foiracfo - — Xk) factors (see Eq. which are to be integrated 
over Xk- Only the term 

r d rj R ej ( rj ) (67) 

remains. Those terms are induced by a vertex in the potential correlation function and may be associated with vertices 
at the positions r*j that carry the momenta (£j,mj). For all the other external momenta, integrations lead to 

dn J dh kfP{k t )R u { n ) Ju (k iri )j £i (hr k ) oc & (r k ) . (68) 

Those terms are easily understood as external lines attached to a vertex at the point r k as they correspond to external 
lines in the potential correlation function. Each external line carries the momenta (ti,mi). One should note that rij 
lines attached to a vertex in the potential correlation function correspond to rij — 1 lines in the temperature correlation 
function plus a vertex which carries an angular momentum. 

Integration over any internal momentum qi will give through -P(<7i)^Dirac(Qi + qj) and Eq. I|66|) : 

£ &(xi,xj)Yt u mYLM(rrj) , (69) 

L,M 

with L = L\= L? k , and M — M| = — ML. These terms are involved as soon as an internal line connects two vertices 
in the potential correlation function. In a diagrammatic representation of the temperature correlation function, they 
may be represented by a double line (one carrying the azimuthal momentum +M, the other carrying —A/) connecting 
the vertices located at xt and Xj . 

Finally, if N is the number of external points and r is the number of vertex connecting n\, ...,n r lines respectively 
in the potential correlation function, the numerical prefactor A a for the temperature correlation function reads 

A a = V ni -l---Vn T -l ^ ' TV > . (7U) 

To summarize, the general temperature ./V-point function may be diagrammatically represented with the following 
rules deduced from the basic diagrammatic rules 

• each of the N harmonic coefficients a£ i Tni is represented cither by an external line or by a vertex with a charge 

• vertices are connected together by internal lines with arbitrary indices attached to them, being (Li, Mi) at one 
end and (L i: —Mi) at the other one, 

• one should then consider all the possible diagram configurations, 
and the following correspondences allow the computation of any diagram 

• each external line with indices (£i,rrii) and attached to a vertex point at Tj contibutes as ^i( r »); 
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• each internal line with indices (Li, Mi) and (Zj, —Mi) connecting two vertices at the points and Tj stands for 



• the weight of a vertex with indices (^,m,) at the point r$ to which rij lines of respective indices 
(4,mi), . . . , (4 4 , mjj are attached is i/ ni i^(r 4 ) (n<+ %7t' 

• if TV is the number of external points and r is the number of vertices the numerical prefactor reads 

, 2 \iV/2 (4tt) n+2t - 2 _ r+5N/2-4^N/2-r-2 
Ui [2^ - ^ ' 7T ' 

• the final value of the diagram is obtained after integration over the radial variables and summation over in the 
internal indices Li and Mi. 

Obviously, if those rules are applied to the diagrams of Figs. ^ [2]and|31 one recovers the expressions l|33|) . I|46|) 
and fBTJl . 

Those rules appear to be useful in the computation of any N-point correlation function in the temperature fluctu- 
ations. They provide compact and simple expressions associated with diagrammatic interpretations. 

Having established the forms of the temperature correlation functions, we pay attention in the following section to 
the behaviors of the bi- and tri-spectra as functions of the multipoles and the configurations. 

IV. SHAPES OF THE CORRELATION FUNCTIONS 

The features of high order correlation functions in the temperature fluctuations are difficult to infer. However, at 
large angular scales, the temperature and potential fluctuations are related in a very simple way because low multipoles 
correspond to very few Fourier modes. Hence, in the Sachs- Wolfe limit (valid up to t ~ 20) the functional forms of 
the temperature and potential correlation functions should be the same. On the contrary, temperature fluctuations 
at small angular scales are induced by numerous Fourier modes of the potential: the initial distribution of potential 
fluctuations is altered by projection effects. This renders the functional form of the temperature correlation functions 
much more intricate. Before we explore general configurations, we first establish the explicit Sachs- Wolfe limits of the 
bi- and tri-spectra and examine specific configurations. 

A. The Sachs- Wolfe limits 

The Sachs- Wolfe limits of the functions Ri , ^ read 

^ SW W = ™W(r-r.) , (71) 
6 

e! w W = ~ cf w , (72) 



whereas the vertex propagator contributes as 



)= I dkeP(k)jl(kn) = ^C™. (73) 



From Eqs. I|71 |l -I|73 |l . the Sachs- Wolfe limit of the reduced bispectrum takes the form 

3/2 



b 



sw 



= -24, 2 g) (^17 + perm.), (74) 



with £1,(2,^3 satisfying the triangular inequalities. As expected, the functional form of the temperature bispectrum 
is the same as the potential bispectrum. 

As the bispectrum scales as ~ C\ at large angular scales, the normalized bispectrum 1)39(1 reduces to a constant in 
the Sachs- Wolfe limit 



24^2 (-) ■ (75) 



Similar calculations can be done for the trispectra. In the line configuration, the reduced trispectrum takes the 
Sachs- Wolfe limit 



«»%± (L) sw = 9^( 27 r) 3 Cf w (Cg* + Cf 2 w ) (Cf 3 w + Cf 4 w ) , (76) 
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normalized bispectrum normalized bispectrum 




5 10 50 100 500 5 10 50 100 500 

1 1 



FIG. 4: Left panel: normalized bispectrum for an equilateral configuration but in a ACDM model (dashed line) and a sCDM 
model (solid line) for V2 — 1. The integrated Sachs Wolfe effect prevents the formation of a plateau at low I in a ACDM model. 
Both curves exhibit acoustic oscillations of roughly the same amplitude although acoutic pics seem to be larger in a ACDM 
model. Right panel: £ 2 (£ + 1) 2 6«£ in a ACDM model (dashed line) and a sCDM model (solid line) for V2 = 1. The amplitude 
of the bispectrum is greater in a ACDM model. Comparison between the two plots shows that the second peak (£ ~ 300) in 
but is only due to the minimum of the Ce's. 



while the L-independent star contribution gives 



star ,t 1 t 2 SW 
I'D- v . 



0,,star /r. \3 /^SW^SWr.SW 



3 other terms) 



(77) 



From Eqs. I|76|l and (|77|l . we recover the functional form of the potential tri-spectra when £±, £2 and L on one hand 
and £3, £4 and L on the other hand satisfy the triangle inequalities. One can note that the scaling of the temperature 
trispectrum at large angular scales depends on the configuration: in the line configuration, the trispectrum scales as 
~ Cf for a given L whereas it scales as ~ Cf in the star configuration. The shape of the temperature tri-spectrum 
then depends on the type of mode couplings that are responsible for the potential tri-spectrum. The normalized 
trispectra defined in Eq. i|62|) reduce in the Sachs- Wolfe limit to 



lineal £2 I T \ 

star 7f\f 2 
Hal* 



9*4 no 



9i/| tar 



(27T) 



(2*)" 



1 



Ci 



(Ce 1 + Ci a ) Cl 13 (Ct a + CeJ {Cf t + Q 4 ) Cl 14 (Ct s + Q 2 ) 



(Ci 1 + Ci 2 ) C L (Ci 3 + Ci A ) [Ct t + Ci 2 ) C L (Ci 3 + C e J 
Cf. Cf, + 3 other terms 



(Cf 1 + Ce 2 ) Cl (Cf 3 + Cf 4 ) 



(78) 
(79) 



B. General case 



In the previous paragraph it has been shown that the structure of the high order correlation functions for the 
temperature field reproduces those of the potential in the Sachs Wolfe limit. In general however this is not the case 
and the temperature correlation functions exhibit acoustic oscillations that depend on the geometric configuration [2(| . 

For example, left panel of Fig 0] shows the normalized bispectrum for an equilateral configuration in a A-cold dark 
matter model (ACDM) and a standard cold dark matter model (sCDM). A Sachs Wolfe plateau is expected in a 
sCDM model while in a ACDM model, the integrated Sachs- Wolfe effect strongly contributes to the bispectrum at 
low £ in such a way that the plateau is not visible any more. For greater t, say £ > 20, we can see the imprints of 
acoustic oscillations on the behaviors of the bispectra which vanish for some values of £. The amplitudes of oscillations 
of the normalized bispectra are comparable. As for the CfS, the period of acoutic oscillations seems to be larger in a 
ACDM model. The right panel of Fig0]shows the bispectrum with a different normalization, i.e. £ 2 (£ + l) 2 bfff. We 
can note the presence of a plateau at the same level for both bispectra in the low £ limit although the amplitude of 
oscillations is greater in a ACDM model. Comparing the two normalizations, we can interpret the secondary peak of 
the normalized bispectrum bui at £ ~ 300 as due to the minimum of the C/s. 

We have similarly computed llne t and stal 't for different configurations. Some of those results are displayed in 
Figs. 03 El and [7] for the star (left panels) and line like mode couplings (center panels). 

Fig-Elcorresponds to the losange configuration £\ = £2 = £3 = £4 = £ whereas Figs.[I)]and[7|respectively correspond 
to the configurations £ x = £2 = L, h = 4 = i + 200 and l x = £ 2 = £, £3 = £ 4 = £ + 300. The different curves 
correspond to different losange shapes: L = £/10, L = £/2, L = £, L = 3£/2 and L = 21. We can see very similar 
behaviors of the trispectra: they all exhibit oscillations with an acoustic peak located at £ ~ 300. The line and star 
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normalized trispectrum normalized trispectrum normalized trispectrum 




FIG. 5: Plot of i\\(L) as a function of I for L = £/W (dashed line), L = 1/2 (dash-dotted line), L = £ (solid red line), L = Z£/2 
(long-dash-dotted green line) and L — 21 (dotted blue line) for v l ^ nc — ;/| tar = 1. The left panel represents the star trispectrum, 
the center panel shows the line trispectrum and the trispectrum from lensing effects is plotted in the right panel. 



normalized trispectrum normalized trispectrum normalized trispectrum 




FIG. 6: Plot of tif 200 e+2oo(L) as a function of I for L — £/10 (dashed line), L = £/2 (dash-dotted line), L = I (solid red line), 
L — 3£/2 (long-dash-dotted green line) and L = 21 (dotted blue line) for i/^ nc = z/jf ar = 1. The left panel represents the star 
trispectrum, the center panel shows the line trispectrum and the trispectrum from lensing effects is plotted in the right panel. 



signals are larger in the low £ range and have roughly the same order of magnitude when 1/3 e ~ j/| tar . In the chosen 
configurations, the trispectra both take positive values at least for I < 500. 

However, the line trispectrum is much larger than the star one for large I as well as for low-L values. A significant 
and specific dependence on the shape of the quadrilateral is clearly visible in these plots and may allow to distinguish 
between the two types of mode couplings. In Figs. El and we can see the L-dependence of the line trispectrum: for 
t < 30, the signal is an increasing function of L while it decreases with L in the range 30 < I < 150. This feature is 
quite different from the star trispectrum whose signal is not a monotonic function of L. 

The chosen configurations exhibit special points: for I ~ 30 for the line like mode coupling and I ~ 80 for the 
star like one, the trispectra take values that do not depend on L. Note that for these points, the value of the 
star trispectrum is very low contrary to the line trispectrum. However, we should stress that we only plotted some 
specific configurations whose features may be quite different from other geometrical configurations. Exploring all the 
configurations remain a daunting task although we expect some geometry to be much more sensitive to primordial 
non-Gaussian signals [T^ . 



C. Weak lensing effects 

At the level of the temperature tri-spectra subtle differences can be observed depending on the nature of the 
potential high order correlation functions. This is then interesting to compare those results with the temperature 
trispectra induce by weak lensing lensing effects. As mentioned in the introduction, this effect is the dominant low 
redshift second order coupling. For parity reason, one expects indeed the trispectrum - rather that the bispectrum 
- to acquire a non negligeable value from lensing effects 0, @, 0> 113 • Lensing effects amount to relate the observed 
temperature contrast, 6T/T, to the primordial one, ST/T, through 

^(7) = ^(7 + *7), (80) 
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normalized trispectrum normalized trispectrum normalized trispectrum 




FIG. 7: Plot of ttfsoo i+3Qo(L) as a function of I for L — £/10 (dashed line), L = 1/2 (dash-dotted line), L = I (solid red line), 
L — 3£/2 (long-dash-dotted green line) and L = 2£ (dotted blue line) for v$ 6 = j/§ tar = 1. The left panel represents the star 
trispectrum, the center panel shows the line trispectrum and the trispectrum from lensing effects is plotted in the right panel. 



where 6j is the lens induced displacement field. Usually, the displacement field is smaller than the angular scale under 
interest so that we may Taylor expand the temperature contrast as 



jr(7 + <*7) = jril) + h ■ V — (7) + o(^7 2 ) 
The first non vanishing contribution to the connected trispectrum is given by 



(81) 



,6T,„ .5T,^ .ST. .ST,. ., ,„ ST . .ST . .. ,„ ST ,„ .ST . ,. . , . . 

(^r(7i)^7(72)^r(73)^7(74)) c = ^ — (71)^(73)) (^ — (72)^(74)} (SjuS^j) + H other terms , (82) 

where the Einstein index summation was used. In the small angle approximation, the four-point function reads 

d 2 l x d 2 £ 2 d 2 £ 3 d 2 £ 4 



\Y-(7i)y"(72)^(73)^(74))c = 



2?r 2tt 2tt 2tt 

d 2 L (p^h + CtJi 
+ (2 «-»3) + (2 ^4) , 



<5Dirac(4+4+^3+4)e ! 



?l "71 +i^2 -72+^3 -73+^4-74 



L 



L 
L 



2 &irac(^l+^2-i) 



where 



C'l 



(2tt)2 



XCMB 7/; 2 fv , | 

dX^T^ 



(83) 



(84) 



Here \ an d ^o(x) are respectively the comoving distance and the comoving angular diameter distance from the 
observer while Pg is the matter power spectrum. The function w is a function of distances and strongly depends on 
the cosmological parameters (see e.g. We define the reduced trispectrum as 



+Cf 



+ct 



C ei £i + C £ J 2 
{Ctji + C ea £ 3 
(CiJx + C e J 4 



L 

L 2 



L13 

^13 



T 2 



and the normalized trispectrum reads 

lcnsT^i^2 ( T \ l cns +^1^2 I T \ 



Ce 2 £2 + Cg 4 £4 
Ci 3 ts + Ci 2 £2 

1 



L 

L 2 



L13 

-^13 
^14 



(Ci 1 + Ci 2 ) Cl (Ce 3 + Ct 4 ) 



(85) 



(86) 



The normalized trispectrum due to lensing effects is plotted in the rig ht p anels of Figs. |3] H3 and [7| for different 
configurations in a sCDM model using the outputs of the CMBFast code [3(J ■ As expected, the contribution of weak 
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FIG. 8: Comparison between the variations of the radial transfer function Rt{r) (dashed line) and £e(r) (solid line) as functions 
of r for different I in a sCDM model. From the upper left panel to the lower right panel: I = 6, I = 50, I = 110, I = 300. 



lensing occurs at small scales and globally increases with I, The signal is enhanced by large values of L whose value 
triggs the positions of the acoustic peaks. 

Weak lensing effects have a specific signature on the temperature trispectrum which is quite different from the 
effects of primordial origin. The value of the trispectrum induced by weak lensing is much larger than the expected 
ones from primordial non-Gaussianities for v l ™ e ~ z/| tar ~ 1. In particular, comparable signals in the range I < 50 
would require 1/3 ~ 10 4 and even 1/3 ~ 10 6 for £ > 200. 

The following sections aim at describing the general behaviors of correlation functions induced by primordial non- 
Gaussian statistics. 



V. PHENOMENOLOGICAL APPROACHES FOR THE RECONSTRUCTION OF THE CORRELATION 

FUNCTIONS 

As we saw in the previous sections, the correlation functions of the coefficients ai m are given in term of three 
functions: the radial transfer function Re(r), the function &(r) and the vertex propagator Q(r,r'). Assuming a given 
power spectrum the propagator Q(r,r') is a known function, peaked at r = r' whereas Rg(r) and £t(r) encode the 
details of the microphysics and have to be computed numerically. The goal of this section is to describe the behaviors 
of these functions in a simple way. As shown in Fig. [8] the variations of £e(r) with r are smoothed compared to 
those of Ry (r) which peaks at the last scattering surface r = r* . Given the shape of the radial transfer function, 
it is then natural to expand it as a combination of a Dirac distribution, that would correspond to the Sachs Wolfe 
limit, and its first derivatives or some equivalently smoothed functions. This expansion seems valid as soon as one 
neglects the late integrated Sachs- Wolfe effect, which is not localized at the last scattering surface. This expansion 
reduces to a description of the temperature anisotropics today through the knowledge of the first radial derivatives 
of the gravitational potential on the last scattering surface: <fii m {r*), 4'£ m (' r *)i etc- This approach should be linked to 
a multipole expansion in the strong coupling limit. The goal of the next subsection is to test the validity of such an 
expansion. 



A. Which expansions for the transfer functions? 



Neglecting the integrated Sachs- Wolfe effect and assuming a nearly instantaneous recombination, the description 
of the baryon-photon plasma dynamics is encoded by the first multipoles Oi(k,ri*) at the last scattering surface. The 
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FIG. 9: Shape of g\ (left panel) and g2 (right panel) as a function of r/r* for I — 10 (solid line) and i — 100 (dashed line). 



transfer function can be approximately written with the following expansion 



(87) 



where, to keep insights of the micro-physics in mind, the monopole and the dipole are very roughly given by expressions 
such that (see |28) for a more accurate phenomenological description of the acoustic pics), 

O (M*) - -^cos(fcc s 77»)e~ fc2/fc o 



i(M*) 



i 



3^ 



sin(/cc s 77*)e 



-k 2 /k 2 L 



(89) 



where the damping effect is schematically taken into account (and kn is the damping scale), 77* is the conformal time 
at which recombination occurs, c s is the sound of speed of the plasma at the last scattering surface. 

In the following we propose to replace the expansion l|87|) . where the coefficients are k dependent, by effective 
^-dependent coefficients, e.g., 



T e (k) = A (£)j e (kr r ) + A 1 (£)j' e (kn) + A 2 (£)j' e '(kn) 



(90) 



This is justified from Eq. I|87|l by the fact that the Bessel functions, and their derivatives, peak at k ~ £/r* . Roughly 
speaking we then expect that Aq(£) « 0o(£/r*,r]*), A%(£) ~ 9\(l/r*, 7?*). This will explain the rough behavior of these 
coefficients. In the following though we will not take this identification into account any more. 

The idea is then, that at least for large enough angular scales, such an expansion could provide a reasonable 
description of the observed anisotropies when only a few terms are taken into account. Such a convergence is ensured 
for Eq. H87[) because the multipole decomposition at the last scattering surface is naturally ordered by a small parameter 
k/f, where f is the recombination rate. Incidentally one can remark that the first term of Eq. I|9U|I would correspond 
to the Sach- Wolfe effect in the limit where Ao is independent on I. In the following though, we treat these coefficients 
on a pure phenomenological footing. We will only assume that the decomposition (|90|l is sensible and provides us 
with a good description of the transfer function. 

Such an expansion implies the following form for the radial transfer function 



Rt (r) = A (£)g (r, n;£) + A 1 (£) 9l (r, r* ; £) + A 2 (£)g 2 (r, r* ; £) 



with 



g n (r,n;£) = 



dk k 2 r 2 



(91) 



(92) 



The c/i function is given in terms of the hypergeometric function 2-F1 and corresponds to a ^-dependent smoothing 
of <5p irac (r — r*). The expressions of g\ and g 2 are given in appendix 1X1 and their shapes are shown for different £ in 
Fig. H Note that g (r, n;t) = 5 Diiac (r - n)n/2^. 

To describe the effect of the expansion (|9Uf) . it is necessary to know the coefficients Ai(£). We then adopt a 
phenomenological approach in which these coefficients are determined from a comparison with the CMBFast code 
results |33. All the following results are obtained in a standard cold dark matter (sCDM) model, where the integrated 
effect may be neglected. 

In principle, one should be capable of extracting the coefficients Aq, A\, A 2 ,... by fitting the numerical radial 
transfer function with a combination of the functions go, g\, g 2 ,... . However, this approach raises some issues. In 
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FIG. 10: Shape of fo (left panel), /i (center panel) and /2 (right panel) as a function of r/n for £ = 10 (solid line) and £ = 20 
(dashed line) (sCDM model). 




FIG. 11: Shapes of the coefficients Aq (left panel), Ai (center panel) and A2 (right panel) as functions of £ in a sCDM model. 
These coefficients have been obtained from a fit of the standard CDM outputs of the CMBFast code with the ansatz Eq. 19311 for 
different choice of n. Solid line: r* = 11670 h" 1 Mpc. Dashed line: n = 11690 h' 1 Mpc. Dotted dashed line: r* = 11700 h _1 
Mpc. 



particular, the determination of the coefficients Ai is extremely sensitive to the choice of r* because the functions gi 
peak at the last scattering surface. On the other hand, the expansion (|90|) induces the following expansion of the 
smoothed functions &(r): 

&(r)=A (£)f (-;l)+A 1 (l)f 1 (-;l) + A 2 (t)f 2 (-;l) + ... , (93) 
where, for a scale-invariant power spectrum, 

fn{ 7j £) = I T Mkr)j(en){kn) ■ (94) 

The functions fi are computed for i = 0,1, 2 in appendix iBl and plotted in Fig.^| Fitting the functions £e(r) with 
a combination of the functions fo, f\ and ji gives an accurate determination of the coefficients Aq, Ai and A 2 . 
The results are shown in Fig. ^2 Comparisons with the expected coefficients A and Ai in Eqs. I|88|l and 189|) for 
typical values of c s and kn are shown in Fig. We can see that the oscillations of the coefficients originate from 
the oscillations of the first multipoles on the last scattering surface. As expected from the Sachs- Wolfe limit, Aq(£) 
reduces to the constant value ~ —1/3 at low £ and A\{1), A 2 (£) — > when £ goes to zero. Note that neglecting the 
early integrated Sachs- Wolfe effect in Eq. (|87|) for a sCDM model do not modify much the expected coefficient Aq at 
low multipole values £. This would not be the case in a ACDM model for instance, where the integrated Sachs- Wolfe 
effect would induce large values of Aq at low I. 



B. Power spectrum 



Given the coefficients A pi the computation of the C/s may be performed analytically. Integrating over the radial 
coordinate r in Eq. (|23)l yields 



C t = -/ ^Ukf , (95) 
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where we set a scale invariant power spectrum. Using the expansion (|90|l leads to 

A\ .„ 3^ + 3^-10 



A\ 



A A 2 



2£(£+l) 6{£~l){£ + 2) z 30£(£ - 2){£ + l){£ + 3) 3£(£+l) 



+A A 



tt(2£-3)11 
1 (2^ + 3)!! 



A X A 



tt(£ 2 +£ + 3)(2£-5)!! 



(2* + 5)!! 



(96) 



where higher order terms were not written explicitly as they are expected to be small corrections to this result. From 
the above expression, we can see that all the terms involving AiAj are subdominant when i + j is odd. This is a quite 
general fact as we will see in the flat sky approximation (see appendix[UJ). For sufficiently large £, only the dominant 
terms contribute and the two-point function becomes 
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A 2 



Al 
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3. 

10 



A A 2 
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(97) 



Using the values of the previously obtained coefficients A, t , Eq. 1)9 7|) appears to be a rough estimate of the power 
spectrum as shown in Fig. 1131 

As we have shown in the previous subsection, the computations of higher order correlation functions involve quite 
cumbersome integrals over hypergeometric functions, which cannot be expressed analytically. To get an approximate 
result, it would seem natural to expand &(r) about r — r» with a Taylor expansion since its behavior is smoothed 
compared to the radial transfer function's. However this series expansion appears not to be valid in a sufficiently 
wide region around r — r». As this approach does not allow to compute high order correlation functions, we should 
use another approximation. For sufficiently large £ (£ > 10 in practice), Fig. 1131 shows that the two-point function is 
well described only by the leading order terms in £. As the large £ approximation corresponds to small angular scales, 
it also means that the sky is assumed to be locally flat so that projection effects reduce to a simple identification. 
In the following we assume that the last scattering surface is properly described by a plan. Hence, the multipolc 
decomposition corresponds to a flat two-dimensional Fourier transform. 

Using the expression of the temperature contrast 
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we get 
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The decomposition of the transfer function in terms of Bessel functions gives 

13 1 



ST 
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(99) 



(100) 



Considering that the functions A p (£) do not vary much compared to the oscillating Bessel functions, we can use 
Eq. (JTSJ) to get 
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1 (V 
kPdri 



with 



r* 



k z u z + k± , 
k z n + £ ■ 7 



(101) 



(102) 
(103) 
(104) 



The x-axis, whose unit vector u z is oriented towards us, supports the line of sight. Vectorial quantities belonging to 
the plan which is orthogonal to the line of sight are indexed by a symbol _L. The vector £ = k±/r* is defined as the 
conjugated variable of the angle 7. 
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FIG. 12: Comparison between the fitted coefficients (soiide fine) Ao(£) (ieft panei) and Ai(£) (right panel) and the roughly 
estimated expected ones (dashed line) from Eqs. I|880 and l|89fl with the typical values c a ~ 1CP 1 , n/r/ t ~ 10 3 and knn ~ 10 3 . 
Eqs. and ® give a rough estimate of the coefficients Ao and and explain their oscillations with I. 



The temperature contrast becomes 
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and the two-point function reads 
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with P(y) = P 



For a scale invariant power spectrum P(y) = 1/y 3 , the following integrals 
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which is in perfect agreement with Eq. <|97|l (see appendix [U] for the correspondence between flat sky and all sky 
formalisms). As we highlighted, only terms AiAj with i + j even do contribute to the two point function. The 
convergence of the series may be checked in Fig. ^| where we used the expansion (|90|l to A4. 

The large £ limit together with the expansion l|90|l provides us with a more tractable formalism to express higher 
order correlation functions. In the following, we explore the approximate behavior of the bispectrum through the very 
first terms of the expansion (|90|l . 



C. Bispectrum 

An accurate description of the bispectrum is a powerful tool towards discriminati ng a mong the inflationary mod- 
els [^IH^ or understanding the primary or secondary nature of non-Gaussianities [U [3 ES Eil • I n this perspective, 
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FIG. 13: Reconstruction of the power spectra from the ansatz l|90|l . The coefficient Ai are obtained from a fit of the standard 
CDM outputs of the CMBFast code. Left panel: effect of the fiat sky approximation on the power spectrum recontruction. 
The solid black line shows the results from CMBFast. The red dashed line shows the results from Eq. Ij96^ whereas the blue 
dotted-dashed line is obtained from the flat sky approximation Eq. 1971 . Those two recontructed curves are in agreement with 
each other for I > 15. Right panel: reconstruction of the power spectrum with different degrees of approximation: expansion to 
first order (contribution of Ao and Ai only) (dotted black line), to second order (long-dashed red line), to third order (dotted 
dashed green line) and to fourth order (short-dashed blue line). The solid black line represents the expected result from the 
CMBFast code. 

we apply our expansion to explore the key features of the temperature three-point function which may be expressed 
as 

!^(7i)^(7 2 )^(73))c= £ J ^-^^ A p^ A p 3 F^^iMe^^e^^e-^S^^+^ + ia) 

(111) 



with 
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perm. 
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From the changes of variables x% — > —Xi, one can see that the non- vanishing terms are those with p\ +P2+P3 even. 
In the special case where p 3 = 0, both p\ and p 2 should be even. 

As the direct computation is uneasy, we decompose the Dirac function fcirac (xi + x 2 + x 3 ) as in Eq. 1|29[) to recast 
the three-point function into 

,ST ST .ST.. .. v 2 fd 2 4 fd 2 £ 2 fd 2 £ 3 /- - -\ _^ 1 . 1i+ ^ +x .t\ u 

(jrhi)jr(l2)jr(j3)}c = ^y^ j J J ferae (4 + 4 + 4 ) e '^''^^'^(HAA) , 

(113) 

with 

Vi A A) = B (<iA,4) + &(t 3 ,tui2) + B(t2,t3,ti) • ( 114 ) 
For a scale invariant power spectrum 
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normalized bispectrum 
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FIG. 14: Normalized bispectrum for an equilateral configuration hut in a sCDM model with 1/3 = 1. The solid black line is 
the expected bispectrum. The red dashed and the green dot-dashed lines are respectively the reconstructed bispectrum to first 
and second order. The blue lines show the behavior of the coefficients Aq(£) (dotted line) and A\(£) (long-dash-dotted line). 
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We defined the following functions (see appendix 
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where K n stands for the modified Bessel function of the second kind of order n. 
The link between the flat sky and all sky formalisms imposes (see appendix Q 
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In the peculiar case where 4 



4 = the functions 4, 4 and 63 simplify to 
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(115) 
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Numerical results are shown in Figs. 1141 1151 and 1161 In Fig. 1141 the normalized bispectrum in l|113|) is compared 
to the expected one in an equilateral configuration. For an equilateral configuration, it appears that the bispectrum 
is properly approximated with a few terms. In particular, one can see that the main features of the bispectrum are 
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normalized bispectrum 
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FIG. 15: Reconstruction of the normalized bispectrum for an isosceles configuration be e 70 as a function of £ in a sCDM model 
with V2 — 1. The solide black line represents the expected bispectrum. The green dotted dashed line and the red dashed line 
represent the reconstructed bispectrum from an expansion to the order Ai and Ai respectively. The blue lines represents the 
coefficients Ao(£) (dotted line) and Ai(£) (long dashed dotted line). 

normalized bispectrum 
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FIG. 16: Normalized bispectrum for an isosceles configuration be e+300 <+300 (left panel) in a sCDM model (solid black line) with 
1/2 = 1. The green dotted dashed line and the red dashed line represent the reconstructed normalized bi-spectrum to order 
Ai(£) and A2(£) respectively. The blue lines represent the coefficients Aq (dotted line) and Ai (long-dash-dotted line). 

encoded by the coefficients A and A\. Recalling Eq. i|115|) . it means that roughly bm oc Aq(£)(Aq(£) + 0.5ylf(£)). 
Hence the zeros of the bispectrum are roughly given by the zeros of the monopole Aq{£) for an equilateral configuration 
(see Fig. ITU. 

Fig El represents the peculiar configuration where one of the I is fixed to a zero of the monopole, £ — 70, while the 
other two are equal. From Eq. 1)115(1 . the bispectrum should reduce to first order to a single term 

&7o«oc4>(ty4i(ty4i(70). (123) 

Hence the zeros of the bispectrum are roughly those of the monopole and the dipole terms. Plots of FigllSlshows that 
this approximation is too crude to estimate the zeros of the bispectrum mainly because this configuration demands 
the dominant term to vanish. However the coarse features are recovered and the sub-dominant terms contribute as a 
global shift that lead to a proper fit of the expected bispectrum. 

We also paid attention to configurations that would correlate the first two acoustic peaks of the C!e's. Fig. 1 161 shows 
the bispectrum 6^+300^+300- The basic features of this bispectrum may be infered by noting that Ao(^+300) ~ —Aq(£) 
and A x {£ + 300) ~ -Ai{£). We then get 

Woo ^+300 *A a (£)[A (£) 2 + . 5A,(£) 2 ] and Wsoo oc -A (£) [A a (£) 2 + h{£)A x [if] , (124) 

where h{€) ~ 0.25 is a slowly varying function. These configurations roughly reproduce the bispectrum corresponding 
to an equilateral configuration. 

Our results provide a simple way towards understanding the behavior of the bispectrum as a function of the 
configuration. Recalling that the coefficients Aq, Ai and Ai represent the first momenta of the gravitational potential 
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on the last scattering surface, we have shown that the coarse features of the bispectrum may be infered from the 
monopole and the dipole behaviors. 



D. Trispectrum 



The analysis of the trispectrum may be performed in the same way as previously for the bispectrum. As the number 
of terms in our phenomenological approach increases with the order of the correlation function, the key features of 
the trispectrum become less obvious to infer. However the sign of the trispectrum may be roughly predicted given 
the behaviors of Ao and A\ . 

The very first terms of the expansion of the trispectrum read 
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For a losange configuration l\ = £ 2 = £ 3 = £4 = £, it may be checked that the signs of trispectra do not change since 
the dominant terms read 
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In a similar way, we could expect that the configurations such as l\ = £ 2 = £, £ 3 = £4 = £' have a constant sign 
for not too large multipoles £ and £'. This seems to be the case in Figs El EH and Fig. 1171 shows the reconstructed 
trispectrum t^f 300 £ +3 qq {L = 5). We can check that the trispectrum does not vanish and that its main features are 
properly reproduced by our expansion in terms of momenta of the gravitational field. 



(125) 



VI. CONCLUSIONS 



We investigated the shapes of high order correlation functions in the CMB anisotropy maps. Formal expressions 
for the high order correlation functions can easily be obtained. We found that, for generic models of primordial 
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normalized trispectrum 
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FIG. 17: Plot of star '7? + 30 (W 30o( L ) as a function of I for L = 5 with i/| tar = 1. The solid line correspond to the CMBFast code 
results whereas the dashed one is the reconstructed trispectrum. 



non-Gaussianities, they can be described with the help of formal diagrams evaluated with computation rules we gave. 
We found that only three kinds of functions with simple diagrammatic interpretations are needed: the radial transfer 
functions Ri{r) represented by vertices, the propagators £i{r) represented by outer lines joining a point to a vertex 
and the vertex propagators, C,L{r;r'), corresponding to internal lines joining two vertices. This formalism provides 
simple expressions of high order connected correlation functions of the observed temperature field for a large class of 
models of primordial non-Gaussianities. 

The development of instabilities at super-horizon scales is intrinsically a nonlinear process that induces specific 
non-Gaussian effects. These types of mode couplings naturally induce a potential bispectrum the amplitude of which 
is characterized by a parameter V2 (e.g. Jnl) of the order unity. However for such a low level of couplings and 
beyond the Sachs Wolfe regime, the calculation of the temperature correlation functions requires a consistent second 
order treatment of the CMB anisotropics through the physics of recombination which is not taken into account in the 
formalism we present. 

More significant sources of non-Gaussianities might come from additional scalar degrees of freedom in the very early 
universe. In such cases, the statistics induced may depend on the details of the scenario and in particular on the types 
of mode couplings during or at the end of the inflationary period. It makes the search of non-Gaussianities of the 
observed CMB maps precious for discriminating between different inflationary scenarii. This is illustrated for instance 
by the qualitative differences we found at the level of the tri-spectra between line- and star-like mode coupling effects. 
Those shapes are furthermore found to be quite different from those induced by weak lensing effects. 

As we stressed it is possible to obtain formally compact expressions for the high order correlation functions. It 
is nonetheless difficult to get insights into the geometrical dependences of these correlation functions in general. 
Projection effects and physics of recombination both contribute to make difficult the transcription from the high 
order correlation functions of the gravitational potential to those of the observed temperature field. Only in the 
Sachs Wolfe regime is this transcription easy. The same formal structure between the correlation functions is indeed 
recovered. To obtain further insights into the behaviors of the high order correlation functions at small angular scale, 
we propose a description of the acoustic oscillations using semi-analytic methods. Ignoring the integrated Sachs- Wolfe 
effect, we found that it was possible to use a gradient-like expansion as suggested by the fact that the radial transfer 
function is peaked in the vicinity of the last scattering surface. This approximation assumes that the temperature 
anisotropics can only be described by the behavior of the potential at or in the vicinity of the last scattering surface. 
In the approach we adopted the transfer physics is encoded by angular scale dependent terms while projection effects 
are treated separately. In this perspective, the small angle approximation is appropriate to handle projection effects 
and leads to more convenient results. We found that the shapes of the temperature power spectrum and bispectrum 
can be reasonably reproduced with the very first terms of such an expansion. 

This perturbative approach can give insights into the behavior of the high order correlation functions. For instance 
it gives a good account of the positions of the zeros of the bispectrum. They indeed have to roughly coincide with those 
of the monopole term at the surface of last scatter. Higher order expansion obviously leads to an increasing accuracy 
in the description of the amplitude of the temperature bispectrum. The trispectra show more complex structures with 
complicated sign patterns but that can also be explained, to a large extent, with the approach we have developed. In 
particular, some key features of the trispectrum, such as its sign in an equilateral configuration, may be accounted 
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for. Nevertheless, the features of the high order temperature spectra might be complementary highlighted in a real 
space treatment of the CMB anisotropies |3l|]. We think that the method we have presented could be extended to 
the properties of the polarization field as well. 
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[33] At best what a complete theory would give is the expected probability distribution function of a computation sketched 
in 

APPENDIX A: FUNCTIONS INVOLVED IN THE EXPANSION OF THE TRANSFER FUNCTION 



Expanding the transfer function as in Eq. 19Uf) induces the calculation of specific integrals in the expression of the 
radial transfer function Ri(r) or in the function £,i{r). In this subsection, we will give the expressions of the function 
gi(r,r*;£) and fi(r/r*;£) (see Eqs. 

To compute these integrals, we will use the following equality (see |2^) for —2 < p <£ + £' + 1 and r < r*: 



dfc k p j/' (kn)je(kr) 



T( e+e '- p+1 ) 



2P+2 r i- P+ i T{ c>-e+ P+2 )r{e + 3/2) 

e+e'-p + i i-£'-v 
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£ + 3/2; 



/ 2 



(Al) 



whereas for r > r* 
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e+e'-p+i 



dk k p ji>(kn,)je(kr) 
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2 ' 2 

On the other hand, the recurrence relation between spherical Bessel functions 

1 



leads to 



i«(*)= £ °$M* 

t'=l-p 



(A2) 



(A3) 



(A4) 



(p) 

where the coefficients a) L satisfy the recurrence relation 



>+*) 



£' + 1 
2£' + 3 



(p) *- (p) 



The integral defining g\ in Eq. I)92|l may be computed with Eqs. (|A1(I - I|A3|I assuming that the limit p - 
analytically continued. The expression for g\ yields for r < 



(A5) 
-2 may be 



Similarly for r > r* 

T(£+l) /r, 
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Here 2F1 stands for the hypergeometric function. In the computation of 52, the integrals in Eq. IjAlfl are not defined 
for p = —2. Hence, we should write g% as 

a2 r 

(A8) 



(A9) 



92(r,n;i) = -j^- / dk je{kr)j e (kn) 



Using Eqs. i|A"Tj) and ifAljl. we get 
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and the function 172 reads 

k (n) for r <r * 
for r > r* 

The expressions of the integrals fi in Eq. I|94fl may be performed with Eqs. ljAl(l - l|A5() and yield for r < 
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In the case r > r» , they transform into 
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APPENDIX B: INTEGRALS INVOLVED IN THE COMPUTATION OF THE BISPECTRUM 



The computation of the bispectrum involves integrals such as 



dx 



{x 2 +^2)(p+3)/2 

I I 

The integrals for p = 0, 1 and 2 may be performed |29j 
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where K n is the modified Bessel function of the second kind of order n and G™^ n ( x 
G-function. The Meijer's G-function is defined by 



ai, a p 
h,...,b q 



a%, a p \ _ 
&i, 



1 



n 3 "Li r fe+ g )n;ur(i-%- S ) 
wtt j c uUn+i r (! - ^ - *) nu+i r ( a J + s ) ' 



(Bl) 

(B2) 
(B3) 
(B4) 

is the Meijer's 
(B5) 
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where the contour C lies between the poles of T(l — dj — s) and the poles of T(bj + s). The last integral (|B4|I may be 
simplified using 

, X (£ 2 + x 2 ) 5 / 2 ~Md^d£ L X + x 2 ) 3 / 2 ■ ( ' 

With Eq. (|B2|) and the differential equations satisfied by the Bessel functions, we obtain 

^ ( / + 6 x 2)5/2 = h W\W\*\) ~ t 2 z 2 K,{t\z\)} . (B7) 
Some other integrals involved in the computation of the bispectrum are of the type 



Results for p = 0, 1 and 2 read 



For the special cases z = 
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APPENDIX C: FLAT SKY FORMALISM 



We define the quantity a(£) by the bidimensional Fourier transform of the temperature contrast 



or conversely 



<■<(<)= I^tt { ' ]v " (C1) 



f (7) = / ^ a(^. (02) 



Inserting the multipole expansion of the temperature contrast in Eq. KTI we may express a(£) as a function of the 
coefficients a£i m > 

/i2 a 
^!W(7)e-^7 (03) 

In the small angle approximation, the spherical harmonics are properly approximated by 

Y im (9,4>)~ J m (19) \[^e im * (C4) 



for I 3> 1 and 6 < 1. Moreover, we can decompose the vectors 7 and i in the small angle approximation 
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The quantity a{jt) then becomes 
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Using the integral representation of the Bessel functions 
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one naturally gets 



which leads to 
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Defining the power spectrum as 



(ag m a}, m ,) = Ct Sa>5 mm >, 



we find, using Eq. (|C9|I that 



(a(/)a(F))=fc irac (l-l ? ) Cfc. 
In an all sky formalism, the reduced bispectrum is defined by 

Using Eqs. HH> and EH, we find the power spectrum in a flat sky formalism 
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The reduced trispectrum is defined by 
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and the correspondence between flat sky and full sky formalisms imposes 
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